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Abstract 

Sum rules provide useful insights into transition strength functions and are often expressed as expectation 
values of an operator. In this letter I demonstrate that non-energy-weighted transition sum rules have strong 
secular dependences on the energy of the initial state. Such non-trivial systematics have consequences: the 
simplification suggested by the generalized Brink-Axel hypothesis, for example, does not hold for most 
cases, though it weakly holds in at least some cases for electric dipole transitions. Furthermore, I show 
the systematics can be understood through spectral distribution theory, calculated via traces of operators 
and of products of operators. Seen through this lens, violation of the generalized Brink-Axel hypothesis 
is unsurprising: one expects sum rules to evolve with excitation energy. Furthermore, to lowest order the 
slope of the secular evolution can be traced to a component of the Hamiltonian being positive (repulsive) or 
negative (attractive). 
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One important way to investigate quantum sys¬ 
tems both experimentally and theoretically are 
through strength functions, 

S{E,,E,) = Y.5{E, - Ef+E,)\{f\d\T)\^.{l) 
f 

which is the probability to transition from a state 
at initial energy Ei to some final state at an en¬ 
ergy Ef = Ei + Ex, via the operator O; |f) and 
I/) are initial and final states, respectively. In par¬ 
ticular I consider strength functions sum rules for 
atomic nuclei for transitions such as electric dipole 
(El), magnetic dipole (Ml), electric quadrupole 
(E2), Gamow-Teller (GT), and others [l|, Q- Such 
transitions not only provide important diagnostics 
of nuclear structure, and thus test theoretical de¬ 
scriptions of nuclei against experiment, but also 
have important impacts in astrophysical physical 
processes such as nucleosynthesis [Sl U, Is], l6|, l7| , neu¬ 
trino transport Q, in the experimental extraction 
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of the density of states i 0, m, El , and so on. 
These applications include initial states excited far 
above the ground state. 

Often one sees the strength function displaying 
either a sharp or a broad peak, which is called a 
resonance, and if most of the strength is in that 
peak, it is a giant resonance [l^. Giant resonances 
can have an intuitive picture: for example, for the 
giant (electric) dipole resonance, or GDR, one en¬ 
visions protons and neutrons collectively oscillating 
against each other [illil- „ _ 

The Brink-Axel hypothesis E113 states that 
if the ground state has a giant electric dipole res¬ 
onance, then the excited states should have giant 
dipole resonances as well; because the GDR is ex¬ 
plained by a collective proton-versus-neutrons os¬ 
cillations, it should not be very sensitive to the de¬ 
tails of the initial state, and so the strength function 
S{Ei,Ex) in ([T]) should independent, or nearly so, 
of Ei 0. While the original Brink-Axel hypothesis 
only concerned the GDR, it later became a simplify- 


^For some historical details, see 
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Figure 1: Non-energy weighted sum rules (here denoted as 
R) as a function of initial energy Ei for computed in the 
sd space. The units for summed Ml strengths (plots (a) and 
(b)) are = nuclear magneton), while those for E2 

(plots (c) and (d)) are e^-fm"^. The (red) dashed line is the 
linear approximation to the secular behavior of R, derived 
from spectral distribution theory, while the (blue) solid line 
is the quadratic approximation. 


ing assumption applied to more general transitions, 
for example Ml and GT. As strength functions off 
excited state are particularly difficult to measure 
experimentally, this hypothesis, if true, would be 
very useful. 

But is the Brink-Axel hypothesis true, especially 
for transitions other than electric dipole? And if it 
is not true, can we do anything about it? 

Despite wide usage and some early experiments 
in support of the Brink-Axel hypothesis [l^, there 
is considerable evidence the Brink-Axel hypothe¬ 
sis fails or is modified for El [ir 
Ml 


21 


E2 21 


mi, 

and Gamow-Teller 
transitions. Nonetheless, as stated in a 
recent Letter [2J , “It is quite common to adopt the 



so-called Brink-Axel hypothesis which states that 
the strength function does not depend on the exci¬ 
tation energy.” On the other hand, a recent ab ini¬ 
tio calculation supported the Brink-Axel hypothesis 
for El transitions from low-lying states [28j and the 
success and consistency of the Oslo method for de¬ 
termining the level density relies upon Ml stren gth s 


following the Brink-Axel hypothesis iEElE 


To simplify the question, I focus on sum rules 
derived from CD, in particular the total strength or 
non-energy-weighted sum rule (NEWSR), 

Rm = j S{E,,E)dE = Y,\{f\d\i)\^ 


= (zidtoK) = mm, (2) 


where for convenience I’ve dehned &0 = R. IfO is 
a non-scalar operator with angular momentum rank 
K and isospin rank I, then R — 

Ok-M, i-fj.OK M,i fj.- (This definition includes a 
sum over charge-changing transitions for 1 = 1, but 
in return i? is a simpler, isoscalar operator; I have 
no reason to believe this qualitatively changes any 
of my conclusions.) Many other sum rules, such 
as the energy-weighted sum rule (EWSR), can also 
be written as expectation values of operators [l^ . 
although here I will only consider the NEWSR. 

If the Brink-Axel hypothesis were true, then 
R{Ei) would be a constant. I investigate the sys- 
tematics of the non-energy-weighted sum rule for 
several different operators O and nuclides as a func¬ 
tion of the initial energy Ei. To test whether or 
not R{Ei) does or does not vary with initial energy 
Ei, I first carry out calculations in a detailed mi¬ 
croscopic model, the configuration-interaction (GI) 
shell model. In the GI shell model, one diag¬ 
onalizes the many-body Hamiltonian in a finite- 
dimensioned, orthonormal basis of Slater deter¬ 
minants, which are antisymmeterized products of 
single-particle wavefunctions, typically expressed in 
an occupation representation [29|. The advantage 
of GI shell model calculations is that one can gen¬ 
erate excited states easily, and for a modest dimen¬ 
sionality one can generate all the eigenstates in the 
model space. 

I use the BIGSTICK GI shell model code [^ . 
which calculates the many-body matrix elements 
Hap = {a\H\P) and then solves H\i) = Ei\i). Greek 
letters {a, (3,...) enote generic basis states, while 
lowercase Latin letters (z,j, .■•) label eigenstates. 
As BIGSTICK computes not only the energies but 
also the wavefunctions, the sum rule R{Ei) is an 
expectation value easy to calculate. 

For this study I use phenomenological spaces 
and interactions, although one could also consider 
ab initio calculations as well; the latter tend to 
have very large dimensions though, making them 
less practical for studying the secular behavior over 
many MeV. Instead I carried out calculations in the 
151/2-0^3/2-0^5/2 or sd shell, u sing a universal sd 
interaction version ‘B’ (USDB) [3l|. I also consider 
the following transition operators: Ml, E2, and 
Gamow-Teller using their standard forms [l, 

I do not use effective charges, I use harmonic oscil¬ 
lator wavefunctions with an oscillator length of 2.5 
fm, I divide sum rules for isovector operators by 
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Figure 2: Non-energy weighted Gamow-Teller sum rules 
(here denoted as il) as a function of initial energy Ei for 
several nuclides in the sd space. The (blue) solid line is 
the quadratic approximation to the secular behavior of R, 
derived from spectral distribution theory. 


3 to roughly average over charge-changing transi- 
tions, and use a quenched value of g a ~ 1; these 
assumptions are scaling factors and do not affect 
my conclusions. 

Fig. [T] shows the NEWSR as a function of ini¬ 
tial energy (relative to the ground state) for the 
nuclide for isoscalar and isovector Ml and E2 
transitions, while Eig. [2] shows the NEWSR for 
Gamow-Teller for several even-even, odd-odd, and 
odd-4. nuclides. (Because I am taking the sum over 
charge-changing transitions and not the difference, 
the Ikeda sum rule will not tumble out of these 
calculations.) I binned the NEWSR into 2 MeV 
bins, but found the size of the fluctuations shown 
by errors bars to be insensitive to the size of the 
bins. Other calculations not shown show qualita¬ 
tively similar results, which can be summarized as: 

• Both the secular (average) behavior of the 
NEWSR and fluctuations thereof show surprisingly 
smooth behavior. 

• As illustrated for Gamow-Teller transitions in 
Eig. [2] (and duplicated but not shown for other op¬ 
erators), the behavior is relatively insensitive to the 
nuclide. 

• Nonetheless, the behavior does depend sharply 
upon the transition: isoscalar E2 falls sharply 
with initial energy, isovector Ml and Gamow-Teller 
grow, and isoscalar Ml has large fluctuations. 

As the Brink-Axel hypothesis predicts NEWSR 
independent of Ei, these results numerically con¬ 
firm the previously mentioned experimental and 


theoretical evidence against Brink-Axel. Given the 
simple yet non-trivial systematics, can one under¬ 
stand these results from basic principles? 

If one is computing R{Ei) in a finite model space, 
such as the Cl shell model, and that model space 
has dimension N, then one can compute the average 
sum rule, that is, the average expectation value, 

1 ^ -I _ 

-^(*|6to|,) = _^(*|i?|,)^(iZ), (3) 

i—1 i 


If we compute the matrix elements of R in some 
orthonormal many-body basis, for example Slater 
determinants in the framework of the Cl shell 
model, the sum is just a trace of the matrix. Be¬ 
cause a trace is invariant under a unitary trans¬ 
formation, we can sum over any convenient set of 
basis states {|q;)}. This invariance under the trace 
is important because the trace can be used as an 
inner product in the space of Hermitian opera¬ 
tors in the framework of spectral distribution the¬ 
ory (SDT), also sometimes called statistical spec¬ 
troscopy QSIlllsilliJlSIli. 

[The notation (R) = {i\R\i) signifies the expec¬ 
tation value being an average over many measure¬ 
ments for the same state. Yet for SDT one aver¬ 
ages the expectation value over all states in a space, 
usually defined by fixed quantum numbers such as 
the number of particles. Practitioners of SDT fre¬ 
quently use the notation for ([3]) H, 37, 3^, 

where m denotes the number of particles and pos¬ 
sibly other quantum numbers, and the trace is im¬ 
plied to be restricted to states with those quantum 
numbers. Because of the unfortunate possibility of 
confusion with the expectation value proper, I in¬ 
troduced a hybrid notation {R) for the average ([3]).] 

To see if the sum rule R{Ei) does indeed have a 
secular dependence upon the initial energy Ei, one 
can take a weighted average, namely, 

^ Y,R{E,)E, = ^Y.^\RH\i)=^). (4) 


N 


Again, since this is a trace, one can compute in 
any convenient basis. Erench proposed [s^ the fol¬ 
lowing inner product between two Hamiltonians, or 
more broadly between two Hamiltonian-like (Her¬ 
mitian and angular momentum scalar) operators: 


= ((i?l - (i?l)) (i?2 - (i?2))) 

= (5) 
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The appeal of this definition of the inner prod¬ 
uct between Hamiltonian-like operators is that, if 
the operators are angular momentum scalars and 
if one works in a finite, spherically symmetric 
shell-model single-particle space, one can calculate 
the traces directly without constructing the matrix 
nisi 113,0. One can sum over states with spec¬ 
ified isospin (while one can take sums over specified 
angular momentum (0 . the resulting formulas are 
significantly more tedious and computationally in¬ 
tensive) or even just on subspaces defined by config¬ 
urations, that is, a fixed number of particles in each 
orbit. In principle one can take higher-order mo¬ 
ments or work with Hamiltonians or particle rank 
higher than two. For this work, however, I use a re¬ 
cent code [0 which reads in only isospin-invariant 
two-body interactions and which calculates at most 
second moments (i.e., the inner product defined 
above) working in spaces with fixed total number 
of valence particles Hvai and total isospin T. 

With the definition of an inner product in the 
space of operators, we can return to the question of 
the invariance of the strength function with initial 
energy. One necessary, but by no means sufficient, 
condition for the invariance of the strength function 
is that the total strength not change with initial 
energy, that is, R{Ei) Ri constant. Such a condition 
implies {RH) {R) (H), but this reduces to the 
inner product {R, H) = {RH) — (R) (H) ~ 0, that 
is, the Hamiltonian H and the operator R must 
be “orthogonal” in a well-defined way. While this 
could happen by accident, in general it will not, as 
we already see in the examples above. 

As it turns out, the above condition corresponds 
to the linear dependence of R{Ei) on Ei. We can go 
to a higher order polynomial description, especially 
if we assume that the state density of the many- 
body Hamiltonian is well-described by a Gaussian 
with centroid E and width a, that is, 


p{E) = N{2Tra- 


2N-1/2 


exp - 


(E-Ef 

2a2 


which is often a good assumption for nuclei 
the language of spectral distribution theory, 

E = W), 




( 6 ) 

. In 

(7) 

( 8 ) 


Let’s further assume that the sum rule R{E) is a 
quadratic polynomial in E: 

R{E) = Ro + R, ^^~^K R2 ~ . (9) 


Then one can easily compute the following averages: 

R = (R) = Ro + R 2 ; (10) 

iR,H)/a = Ri. (11) 


One could add an additional constraint by by higher 
moments, for example {RIE^). While such higher 
moments are calculable [0 . the formula are cum¬ 
bersome and prone to slow evaluation; further¬ 
more experience in unpublished work suggests even 
higher moments have difficulty in describing the 
tails of distributions (This is understandable; 
the traces are just averages, after all, and dom¬ 
inated by the density of states in the middle of 
the spectrum.) Instead, I use the sum rule at the 
ground state energy Egg, which is often accessible: 

R{Egg) = 


These equations (I10I11I12I) can be easily solved. 

Fig [1] shows both linear (red dashed lines) 
and quadratic (blue solid lines) approximations to 
R{Ei). Although the linear approximation demon¬ 
strates a secular dependence on Ei, in general the 
quadratic does better in describing the secular evo¬ 
lution of the sum rule. Fig. [2] shows only the 
quadratic approximation. 

Now, as illustrated in the figures, while one has 
smooth secular behavior, there are nontrivial fluc¬ 
tuations about the secular trends. The fluctuations 
are insensitive to the size of the energy bins. Al¬ 
though the fluctuations about the smooth secular 
behavior are not easily written in terms of traces, 
one might be able to derive the fluctuations from 
random matrix theory; but this will have to be left 
to future work. 

The original Brink-Axel hypothesis described El 
strength functions. To explore them, I use a space 
with opposite parity orbits, the 0 pi/ 2 - 0 p 3 / 2 -lsi/ 2 - 
Od ^/2 or p-sd ^/2 space, chosen so I could fully di¬ 
agonalize for some nontrivial cases. For an inter¬ 
action I use the Cohen-Kurath (CK) matrix ele¬ 
ments in the Op shellim , the older USD interac¬ 
tion ^ in the 0d5/2-lsi/2 space, and the Millener- 
Kurath (MK) p-sd cross-shell matrix elements (i^. 
Within the p and sd spaces I use the original spac¬ 
ing of the single-particle energies for the CK and 
USD interactions, respectively, but then shift the 
sd single-particle energies up or down relative to 
the p-shell single particle energies to get the first 3“ 
state at approximately 6.1 MeV above the ground 
state. The rest of the spectrum, in particular the 
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Figure 3: Non-energy weighted isovector El sum rules (here 
denoted as R) as a function of initial energy Ei for several 
nuclides in the p-sd^i 2 space. The (blue) solid line is the 
quadratic approximation to the secular behavior of i?, de¬ 
rived from spectral distribution theory. 

first excited 0+ state, is not very good, but the idea 
is to have a non-trivial model, not exact reproduc¬ 
tion of the spectrum. Because this space does not 
allow for exact center-of-mass projection I restrict 
myself to isovector El transitions. The resulting 
NEWSRs are shown in Fig. [3l illustrating only a 
weak violation of Brink-Axel. 

Although the quadratic approximation captures 
the general trends, the secular behavior for R{E) is 
not as smooth. This may be because of the model 
space. The density of states for these nuclides, for 
example, are not as Gaussian-like as for the sd- 
shell examples shown; the beryllium and boron nu¬ 
clides have large third moments, while the the neon, 
sodium, and aluminum nuclides have larger fourth 
moments (“fat tails”) than Gaussians. 

Nonetheless, not only do we have evidence that 
the generalized Brink-Axel hypothesis is not fol¬ 
lowed, we can understand why. Previous work has 
suggested specific reasons for breaking the Brink- 
Axel hypothesis: changes in deformation as one 
goes up in e nerg y explains the increase in width 
for the GDR [20|, while a decrease in spatial sym¬ 
metry/increase in SU(4) symmetry explains the in¬ 
crease of strength in Gamow-Teller sum rules [^ . 
Spectral distribution theory provides a more gen¬ 
eral understanding. By establishing a vector space 
for Hamiltonians such that 

= (13) 


the inner product defined by SDT yields {H, R^) = 
Ca (up to some easily-defined normalization). Here 
is the key point: if Ca < 0, that is, attractive, 
one expects a negative slope to R{E') and more 
strength for low-lying initial states. This is seen 
in Fig ([T])(c), where the operator R ^ Q ■ Q, the 
quadrupole-quadrupole interaction. If, on the other 
hand, > 0, that is repulsive, as for (ctt)^ as in 
Fig. ([2]), low-lying states have less total strength. 
Only if Cct ~ 0 could the Brink-Axel hypothesis be 
true, at least at the lowest level. Of course, the lin¬ 
ear approximation is not always sufficient to fully 
describe the secular behavior; for many cases one 
needs at lea st q uadratic and possibly even higher- 
order terms |40|. 

In summary, I have numerically demonstrated 
that that the non-energy weighted sum rule for 
transition operators applied to several sample nu¬ 
clides evolves with the energy of the initial state- 
weakly for the case of isovector EI, and more 
strongly for other operators-and furthermore that 
such variation is expected from spectral distribution 
theory. In particular, one can predict qualitatively 
whether a sum rule will grow or shrink in magni¬ 
tude with initial energy, depending if part of the 
Hamiltonian (that part proportional to the opera¬ 
tor R for the sum rule, that is, the square of the 
transition operator for the NEWSR) is attractive 
or repulsive. In many cases one needs higher mo¬ 
ments for accurate quantitative predictions , but it 
should be clear now that one should only invoke 
Brink-Axel with caution. 

This material is based upon work supported by 
the U.S. Department of Energy, Office of Science, 
Office of Nuclear Physics, under Award Number 
DE-FG02-96ER40985. 
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